FERMI NATIONAL ACCELERATOR LABORATORY 



FERMILAB-TM-2347-E 
TEVEWWG/top 2006/01 
CDF Note 8162 
D0 Note 5064 

18-March-2006 



Combination of CDF and D0 Results 
on the Mass of the Top Quark 



The Tevatron Electroweak Working Group ^ 
for the CDF and D0 Cohaborations 



Abstract 

We summarize the top-quark mass measurements from the CDF and D0 experiments at 
Fermilab. We combine published Run-I (1992-1996) measurements with the most recent 
prehminary Run-11 (2001-present) measurements using up to 750 pb~^ of data. Taking 
correlated uncertainties properly into account the resulting preliminary world average 
mass of the top quark is Mt = 172.5 it 1.3(stat) it 1.9(syst) GeV/c^, which corresponds to 
a total uncertainty of 2.3 GeV/c'^. The top-quark mass is now known with a precision of 
1.3% — a 20% improvement relative to the previous combination. 
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1 Introduction 



The experiments CDF and D0, taking data at the Tevatron proton-antiproton coUider located 
at the Fermi National Accelerator Laboratory, have made several direct experimental measure- 
ments of the top-quark pole mass, M^. The pioneering measurements were based on about 
100 pb"^ of Run-I (1992-1996) data |I[|^ElllliniEllZlEllllinilIIl[I2|and include results from 
the tt qq'bqq'b (all-j), the ti — > iuqq'bb (1+j), and the ti — > £~^h'bi~Vb (di-1) decay channels^. 
Results using approximately 350 pb~^ of Run-II (2001-present) data have been recently pub- 
lished in the l-|-j and di-1 channels [Til E] • The Run-II measurements summarized here are 
the most recent preliminary results in the l-|-j and di-1 channels using 370 - 750 pb"^ of data 
and improved analysis techniques [ini HZl UHl CHI • 

This note reports the world average top-quark mass obtained by combining five published 
Run-I measurements [21 El El 13 QHl ^] with the most recent preliminary Run-II measurements 
from CDF (THl (HI and D0 [THllinj. The combination takes into account the statistical and 
systematic uncertainties and their correlations using the method of references |20| 1^ and 
supersedes previous combinations [2211231 • The most precise individual measurements of Mt are 
now the preliminary measurements in the l-|-j channel from Run II. These are 173.4±2.8 GeV/c^ 
(CDF, jini) and 170.6 ±4.6 GeV/c^ (D0, |TH|). These have weights in the new Mt combination 
of 58% and 21%, respectively. 

The input measurements and error categories used in the combination are detailed in Sec- 
tion 121 and El respectively. The correlations used in the combination are discussed in Section |31 
and the resulting world average top-quark mass is given in Section El A summary and outlook 
are presented in Section El 



2 Input Measurements 

For this combination nine measurements of Mt are used, five published Run-I results and 
four preliminary Run-II results. In general, the Run-I measurements all have relatively large 
statistical uncertainties and their systematic uncertainty is dominated by the total jet energy 
scale (JES) uncertainty. In Run-II both CDF and D0 take advantage of the larger tt samples 
available and employ new analysis techniques to reduce both these uncertainties. In particular 
the JES is constrained using an in-situ calibration based on the invariant mass of W qq' 
decays in the 1+j channel. The Run-II D0 analysis in the l-|-j channel constrains the response of 
light-quark jets using the in-situ W — >■ qq' decays. The JES determined in this manner is then 
propagated to their measurement in the di-1 channel. Residual JES uncertainties associated 

^Here £ = e or /i. Decay channels with exphcit tau lepton identification are presently under study and are 
not yet used for measurements of the top-quark mass. 
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with 7]— and pT-dependencies as well as uncertainties specific to the response of 6-jets are 
treated separately. Similarly, the Run-II CDF analysis in the 1+j channel also constrains the 
JES using the in-situ W qq' decays. The JES is also determined using "external" calibration 
samples as was done for the Run-I measurements. This external JES is applied as an additional 
constraint in this CDF analysis to improve further the total JES uncertainty. Small residual 
JES uncertainties arising from 77— and p^-dependencies and the modeling of 6-jets are included 
in separate error categories. The Run-II CDF measurement in the di-1 channel uses only the 
externally determined JES, some parts of which are correlated with the externally determined 
Run-I JES as noted below. 

The CDF Run-II measurement in the l-|-j channel requires special treatment in order to ac- 
count more accurately for the correlations of the JES uncertainties since the fit uses information 
from both the in-situ and external JES calibrations. In the combination we treat this one mea- 
surement as two separate inputs - one which includes only the in-situ JES calibration, (l+j)i, 
and a second which includes only the JES as determined from the external calibration samples, 
(l-|-j)e. We correlate the JES related error categories as described below while taking the rest of 
the error categories to be 100% correlated between these two inputs. The combination of just 
these two inputs using these correlations yields the identical central value, statistical, JES, and 
total systematic uncertainty as the measurement reported in reference jTH]. The correlations 
between these two inputs and the rest of the inputs are as described in Section HJ 

The inputs used in the combination are summarized in Table Q with their uncertainties 
sub-divided into the categories described in the next Section. 



3 Error Categories 

We employ the same error categories as used for the previous world average ^23j. They have 
evolved to include a detailed breakdown of the various sources of uncertainty and aim to lump 
together sources of systematic uncertainty that share the same or similar origin. For example, 
the "Signal" category discussed below includes the uncertainties from ISR, FSR, and PDF - 
all of which affect the modeling of the tt signal. Additional categories have been added in 
order to accommodate specific types of correlations. For example, the jet energy scale (JES) 
uncertainty is sub-divided into several components in order to more accurately accommodate 
our best estimate of the relevant correlations. Each error category is discussed below. 

Statistical: The statistical uncertainty associated with the Mt determination. 

iJES: That part of the JES uncertainty which originates from in-situ calibration procedures 
and is uncorrelated among the measurements. In the combination reported here it cor- 
responds to the statistical uncertainty associated with the JES determination using the 
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Run-I published 


Run-Il preliminary 






CDF 




D0 






CDF 




D0 






all-j 


1+j 


di-1 


1+j 


di-1 




(l+j)e 


di-1 


1+j 


di-1 


Lumi (pb ^) 


110 


105 


110 


125 


125 


680 


680 


750 


370 


370 


Result 


186.0 


176.1 


167.4 


180.1 168.4 


173.4 


173.4 


164.5 


170.6 176.6 


iJES 


0.0 


0.0 


0.0 


0.0 


0.0 


2.2 


0.0 


0.0 


0.0 


0.0 


aJES 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.8 


1.2 


bJES 


0.6 


0.6 


0.8 


0.7 


0.7 


0.6 


0.6 


0.5 


0.7 


1.2 


cJES 


3.0 


2.7 


2.6 


2.0 


2.0 


0.0 


2.0 


2.2 


0.0 


0.0 


dJES 


0.3 


0.7 


0.6 


0.0 


0.0 


0.0 


0.6 


0.8 


3.4 


3.1 


rJES 


4.0 


3.4 


2.7 


2.5 


1.1 


0.0 


2.2 


1.1 


0.0 


0.0 


Signal 


1.8 


2.6 


2.8 


1.1 


1.8 


0.7 


0.7 


0.8 


0.5 


1.2 


BG 


1.7 


1.3 


0.3 


1.0 


1.1 


0.5 


0.5 


1.2 


0.4 


0.2 


Fit 


0.6 


0.0 


0.7 


0.6 


1.1 


0.8 


0.8 


0.8 


0.5 


0.7 


MC 


0.8 


0.1 


0.6 


0.0 


0.0 


0.3 


0.3 


0.5 


0.0 


0.0 


UN/MI 


0.0 


0.0 


0.0 


1.3 


1.3 


0.0 


0.0 


0.0 


0.0 


0.0 


Syst. 


5.7 


5.3 


4.9 


3.9 


3.6 


2.6 


3.3 


3.1 


3.7 


3.8 


Stat. 


10.0 


5.1 


10.3 


3.6 


12.3 


1.7 


1.7 


4.5 


2.8 


11.2 


Total 


11.5 


7.3 


11.4 


5.3 


12.8 


3.1 


3.7 


5.5 


4.6 


11.8 



Table 1: Summary of the measurements used to determine the world average M^. As described 
in the text, the CDF Run-11 measurement in the lepton+jets channel is treated as two inputs 
in order to more accurately account for the correlations in the jet energy scale uncertainties. 
All numbers are in GeV/c^. The error categories and their correlations are described in the 
text. The total systematic uncertainty and the total uncertainty are obtained by adding the 
relevant contributions in quadrature. 

W qq' invariant mass in the CDF Run-Il 1+j measurement. Residual JES uncertain- 
ties, which arise from effects not considered in the in-situ calibration, are included in 
other categories. 

aJES: That part of the JES uncertainty which originates from differences in detector e/h 
response between 6-jets and light-quark jets. It is specific to the D0 Run-11 measurements 
and is taken to be uncorrelated with the D0 Run-I and CDF measurements. 

bJES: That part of the JES uncertainty which originates from uncertainties specific to the 
modehng of 6-jets and which is correlated across all measurements. For both CDF and 
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D0 this includes uncertainties arising from variations in the semi-leptonic branching 
fraction, 6-fragmentation modehng, and differences in the color flow between 6-jets and 
light-quark jets. These were determined from Run-II studies but back-propagated to the 
Run-I measurements, whose rJES uncertainties (see below) were then corrected in order 
to keep the total JES uncertainty constant. 

cJES: That part of the JES uncertainty which originates from modeling uncertainties corre- 
lated across all measurements. Specifically it includes the modeling uncertainties associ- 
ated with light-quark fragmentation and out-of-cone corrections. 

dJES: That part of the JES uncertainty which originates from limitations in the calibra- 
tion data samples used and which is correlated between measurements within the same 
data-taking period (ie. Run 1 or Run II) but not between experiments. For CDF this 
corresponds to uncertainties associated with the ry-dependent JES corrections which are 
estimated using di-jet data events. For D0 Run-II this corresponds to uncertainties as- 
sociated with the light-quark response as determined using the W — * qq' invariant mass 
in the 1+j channel and propagated to the di-1 channel. The residual ^^-dcpcndent and 
PT-depcndcnt uncertainties for the D0 Run-II measurements are also included here since 
they are constrained using Run-II 7+jet data samples. 

rJES: The remaining part of the JES uncertainty which is correlated between all measurements 
of the same experiment independent of data-taking period, but is uncorrelated between 
experiments. This is dominated by uncertainties in the calorimeter response to light- 
quark jets. For CDF this also includes small uncertainties associated with the multiple 
interaction and underlying event corrections. 

Signal: The systematic uncertainty arising from uncertainties in the modeling of the ti signal 
which is correlated across all measurements. This includes uncertainties from variations 
in the ISR, FSR, and PDF descriptions used to generate the tt Monte Carlo samples 
that calibrate each method. It also includes small uncertainties associated with biases 
associated with the identification of 6-jets. 

Background: The systematic uncertainty arising from uncertainties in modehng the domi- 
nant background sources and correlated across all measurements in the same channel. 
These include uncertainties on the background composition and shape. In particular 
uncertainties associated with the modeling of the QCD multi-jet background (all-j and 
l-|-j), uncertainties associated with the modeling of the Drell-Yan background (di-1), and 
uncertainties associated with variations of the fragmentation scale used to model W-|-jets 
background (all channels) are included. 

Fit: The systematic uncertainty arising from any source specific to a particular fit method, 
including the finite Monte Carlo statistics available to cahbrate each method. 

Monte Carlo: The systematic uncertainty associated with variations of the physics model 
used to calibrate the fit methods and correlated across all measurements. For CDF it 
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includes variations observed when substituting PYTHIA I2S1 12n| (Run I and Run II) 
or ISAJET EH (Run I) for HERWIG |2Hl 121 when modehng the ti signal. Similar 
variations are included for the D0 Run-I measurements. The D0 Run-II measurements 
use ALPGEN jSUj to model the ti signal and the variations considered are included in 
the Signal category above. 

UN/MI: This is specific to D0 and includes the uncertainty arising from uranium noise in the 
D0 calorimeter and from the multiple interaction corrections to the JES. For D0 Run-I 
these uncertainties were sizable, while for Run-II owing to the shorter integration time 
and in-situ JES determination, these uncertainties are negligible. 



These categories represent the current preliminary understanding of the various sources of 
uncertainty and their correlations. We expect these to evolve as we continue to probe each 
method's sensitivity to the various systematic sources with ever improving precision. Varia- 
tions in the assignment of uncertainties to the error categories, in the back-propagation of the 
bJES uncertainties to Run-I measurements, in the approximations made to symmetrize the 
uncertainties used in the combination, and in the assumed magnitude of the correlations all 
negligibly effect (-C lOOMeV/c^) the combined Mt and total uncertainty. 



4 Correlations 



The following correlations are used when making the combination: 



• The uncertainties in the Statistical, Fit, and iJES categories are taken to be uncorrected 
among the measurements. 

• The uncertainties in the aJES and dJES categories are taken to be 100% correlated among 
all Run-I and all Run-II measurements on the same experiment, but uncorrelated between 
Run I and Run II and uncorrelated between the experiments. 

• The uncertainties in the rJES and UN/MI categories are taken to be 100% correlated 
among all measurements on the same experiment. 

• The uncertainties in the Background category are taken to be 100% correlated among all 
measurements in the same channel. 

• The uncertainties in the bJES, cJES, Signal, and Generator categories are taken to be 
100% correlated among all measurements. 



Using the inputs from Tableland the correlations specified here, the resulting matrix of total 
correlation co-efficients is given in Table El 
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Run-I published 


Run-II preliminary 


CDF 
1+j di-1 all-j 


D0 
1+j di-1 


CDF 

(1+j), (l+j)e di-1 


D0 
1+j di-1 


CDF-I 1+j 
CDF-I di-1 
CDF-I all-j 


1.00 

0.29 1.00 

0.32 0.19 1.00 








D0-I 1+j 
D0-I di-1 


0.26 0.15 0.14 
0.11 0.08 0.07 


1.00 

0.16 1.00 






CDF-II (1+j), 
CDF-II (1+j), 
CDF-II di-1 


0.13 0.07 0.05 
0.57 0.32 0.39 
0.30 0.19 0.21 


0.10 0.04 
0.29 0.12 
0.19 0.11 


1.00 

0.41 1.00 

0.06 0.40 1.00 




D0-II 1+j 
D0-II di-1 


0.07 0.04 0.03 
0.04 0.03 0.02 


0.06 0.02 
0.03 0.02 


0.07 0.06 0.03 
0.04 0.04 0.03 


1.00 

0.24 1.00 



Table 2: The resulting matrix of total correlation coefficients used to determined the world 
average top quark mass. 

The measurements are combined using a program implementing a numerical minimiza- 
tion as well as the analytic BLUE method [201 1^- The two methods used are mathematically 
equivalent, and are also equivalent to the method used in an older combination and give 
identical results for the combination. In addition, the BLUE method yields the decomposition 
of the error on the average in terms of the error categories specified for the input measure- 
ments pij . 



5 Results 



The combined value for the top-quark mass is: 

Mt = 172.5 ± 2.3 GeV/c^ , (1) 

with a of 8.1 for 8 degrees of freedom, which corresponds to a probability of 42% indicating 
good agreement among all the input measurements. The total uncertainty can be sub-divided 
into the contributions from the various error categories as: Statistical (±1.3), total JES (±1.6), 
Signal (±0.7), Background (±0.5), Fit (±0.4), Monte Carlo (±0.2), and UN/MI (±0.2), for a 
total Systematic (±1.9), where all numbers are in units of GeV/c^. The pull and weight for each 
of the inputs are listed in Table El The input measurements and the resulting world average 
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Run-I published 


Run-II preliminary 


CDF 

1+j di-1 all-j 


D0 
1+j di-1 


CDF 

(1+j), (l+j)e di-1 


D0 
1+j di-1 


Pull 

Weight [%] 


+0.52 -0.46 +1.20 
-1.6 +0.1 +0.2 


+ 1.59 -0.33 
+ 10.6 +1.1 


+0.43 +0.30 -1.61 
+45.0 +11.8 +10.8 


-0.48 +0.35 
+21.0 +1.0 



Table 3: The pull and weight for each of the inputs used to determine the world average mass 
of the top quark. See Reference [201 ^ discussion of negative weights. 

mass of the top quark are summarized in Figure H 

In the previous combination, for the CDF Run-II 1+j measurement using 320 pb^^of data, 
the in-situ and external JES calibrated inputs each carried approximately the same weight. In 
the combination reported here, the weight of the CDF Run-II 1+j input using the in-situ JES 
calibration carries three times the weight of its counterpart using the external JES calibration. 
This trend is expected to continue since the in-situ JES uncertainty is expected to improve as 
the statistics of the W qq' sample increase with larger data sets. In contrast the uncertainty 
on the external JES calibration already has large contributions from modeling uncertainties 
which may not be reduced with larger data sets. 

The weight of the CDF Run-I 1+j measurement is negative. In general, this situation can 
occur if the correlation between two measurements is larger than the ratio of their total uncer- 
tainties. This is indeed the case here, where the CDF Run-I 1+j measurement is 57% correlated 
with the externally calibrated CDF Run-II 1+j measurement, but with a total uncertainty that 
is about twice as large. In these instances the less precise measurement will usually acquire a 
negative weight. While a weight of zero means that a particular input is effectively ignored in 
the combination, a negative weight means that it affects the resulting central value and helps 
reduce the total uncertainty. See reference [20] for further discussion of negative weights. 

Although the from the combination of all measurements indicates that there is good 
agreement among them, and no input has an anomalously large pull, it is still interesting to 
also fit for the top-quark mass in the all-j, 1+j, and di-1 channels separately. We use the same 
methodology, inputs, error categories, and correlations as described above, but fit for the three 
physical observables, Mf^~'\ M]^\ and Mf^~^. The results of this combination are shown in 
Table [Hand have of 3.6 for 6 degrees of freedom, which corresponds to a probability of 73%. 
These results differ from a naive combination, where only the measurements in a given channel 
contribute to the determination in that channel, since the combination here fully accounts 
for all correlations, including those which cross-correlate the different channels. Using the 
results of Tabled] we calculate the chi-squared consistency between any two channels, including 
all correlations, as x^{.dil — Ij) = 3.0, ~ allj) = 1.3, and x^(c^^^J ~ "^^O — 3-0. These 
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Parameter 


Value (GeV/c^) 


Correlations 




185.4 ± 10.7 
173.4 ± 2.4 
166.1 ± 4.3 


1.00 

0.20 1.00 

0.11 0.32 1.00 



Table 4: Summary of the combination of the nine measurements by CDF and D0 in terms of 
three physical quantities, the mass of the top quark in the all-jets, lepton+jets, and di-lepton 
channel. 

correspond to chi-squared probabilities of 8%, 25%, and 8%, respectively, and indicate that the 
determinations of from the three channels are consistent with one another. 



6 Summary 

A preliminary combination of measurements of the mass of the top quark from the Teva- 

tron experiments CDF and D0 is presented. The combination includes five published Run-I 
measurements and four preliminary Run-II measurements. Taking into account the statistical 
and systematic uncertainties and their correlations, the preliminary world-average result is: 
Mt = 172.5 ± 2.3 GeV/c^. 

The mass of the top quark is now known with an accuracy of 1.3%, limited by the systematic 
uncertainties, which are dominated by the jet energy scale uncertainty. This systematic is 
expected to improve as larger data sets are collected since new analysis techniques constrain 
the jet energy scale using in-situ W —>■ qq' decays. It can be reasonably expected that with the 
full Run-II data set the top-quark mass could be known to better than 1%. To reach this level of 
precision further work is required to determine more accurately the various correlations present, 
and to understand more precisely the 6-jct modeling, Signal, and Background uncertainties 
which may limit the sensitivity at larger data sets. Limitations of the Monte Carlo generators 
used to calibrate each fit method may also become important as the precision reaches the 1% 
level and will warrant further study in the future. 
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Mass of the Top Quark (*Preliminary) 
Measurement 



CDF-I di-l 
D0-I di-l ■ 
CDF-II di-l* 
D0-II di-r 
CDF-I l+j 
D0-I l+j 
CDF-II l+j* 
D0-II l+j* 
CDF-I all-j 



Tevatron Run- 



150 



M 



M^^p [GeV/cl 
167.4 ± 11.4 
1 68.4 ± 12.8 
164.5± 5.5 
176.6 ± 11.8 
176.1 ± 7.3 
180.1 ± 5.3 
173.4± 2.8 
170.6± 4.6 
- 186.0 ±11.5 



XVdof = 8.1/8 



172.5± 2.3 



too 



170 

[GeV/c^] 



190 



Figure 1: A summary of the input measurements and resulting world average mass of the top 
quark. 
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